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SUMMARY 


Synthesis  procedures  for  nonlinear  stochastic  feedback  control  of  gas 
turbine  engines  were  developed  and  evaluated.  Modern  estimation  and  control 
procedures  based  upon  separating  the  stochastic  and  deterministic  aspects  of 
the  control  problem  were  employed.  The  resulting  closed-loop  control  con- 
sists of  (l)  nonlinear  deterministic  feedback  control  logic  designed  using 
piecewise-linear/piecewise-optimal  techniques,  and  (2)  an  estimator  designed 
using  nonlinear  filtering  logic.  Engine  variables  estimated  from  noise- 
corrupted  measurements  are  fed  back  through  the  deterministic  control  logic 
to  generate  commanded  inputs  to  the  engine.  Mode-switching  logic  was  devel- 
oped to  provide  smooth  transition  between  small-signal  regulation  and  large- 
signal  transient  modes  of  estimator /controller  operation. 

The  combined  deterministic  controller/stochastic  estimator  was  evaluated 
by  application  to  a nonlinear  digital  dynamic  computer  simulation  of  the 
FlOO/FUOl  turbofan  engine.  Stability-in-the-large  as  well  as  stability-in- 
the-small  at  design  and  off-design  operating  conditions  was  established. 

Engine  performance  throughout  the  idle  to  military  sea-level  static  operating 
regime  (9  to  100  percent  thrust)  was  investigated.  Thrust  response  of  the 
nonlinear  F100/FU01  engine  simulation  was  evaluated  for  (l)  nominal  and 
degraded  engine  models,  (2)  nominal  and  off-nominal  noise  statistics  and 
(3)  small-signal  as  well  as  large-signal  time  responses.  Stochastic  engine 
performance  (feedback  of  estimated  measurable  and  unmeasurable  engine  vari- 
ables through  the  deterministic  control  logic)  was  compared  with  deterministic 
engine  performance  (feedback  of  actual  engine  variables  through  the  control 
logic).  In  addition,  computational  requirements  of  the  nonlinear  controller/ 
estimator  were  delineated. 
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The  Naval  Air  Development  Center  study  of  Ref.  1 — dealing  with  dynamic 
response  requirements  for  V/STOL  lift/cruise  engines  used  to  provide  height 
and  moment  control  — highlights  a problem  of  increasing  concern  to  the 
Naval  aviation  community.  This  study  points  out  the  inability  of  present  day 
gas  turbine  lift/cruise  engines  — when  used  directly  for  differential  moment 
and  force  generation  — to  meet  acceleration  and  height  control  requirements 
of  V/STOL  Flying  Qualities  Specification  MIL-F-83300.  Current  control  of 
jet-lift  V/STOL  aircraft  at  hover  and  low  speeds  is  therefore  accomplished 
by  employing  an  auxiliary  reaction  control  system  which  diverts  a portion  of 
main  engine  airflow  to  the  aircraft  outer  extremities  to  generate  the  appro- 
priate flight  control  moments  and  forces.  Hover  and  low  speed  control  is 
accomplished  by  modulation  of  airflow  nozzles  which  can  be  rapidly  opened  or 
closed.  This  mode  of  auxiliary  control  results  in  increased  aircraft  weight 
and  complexity  as  well  as  decreased  internal  volume  available  for  fuel.  If 
the  main  engines  could  be  used  to  provide  the  required  rapid  acceleration 
directly,  considerable  weight  and  volume  savings  would  be  realized.  These 
considerations  will  assume  even  greater  importance  in  the  future  as  the  Navy 
progresses  with  development  of  V/STOL  A and  B aircraft  concepts. 

Frevious  UTRC ' programs  with  ONR  have  demonstrated  the  ability  of  multi- 
variable  feedback  control  logic  — using  feedback  of  actual,  not  estimated, 
variables  — to  provide  significantly  more  rapid  propulsion  system  dynamic 
response  than  that  associated  with  conventional  control  logic  (Refs.  2 and  3). 
The  Ref.  2 report  documents  this  improvement  in  dynamic  response  with  partic- 
ular application  to  the  Navy  XFV-12A  V/STOL  aircraft.  Results  obtained  in  a 
UTRC  study  supported  by  the  National  Science  Foundation  indicate  that  use  of 
adaptive  multivariable  control  logic  for  static  performance  optimization  leads 
to  reduced  thrust  specific  fuel  consumption  (Ref.  h).  However,  the  UTRC 
multivariable  control  logic  of  Refs.  2,  3 and  h employs  feedback  of  actual, 
not  estimated,  engine  variables.  In  practice,  key  engine  variables  are  either 
unavailable  for  feedback  or  are  measured  through  noisy,  imperfect  sensors. 

A practical  design  approach  must  therefore  achieve  rapid  propulsion  system 
response  while  accounting  for  (l)  unavailability  of  key  output  variables, 

(2)  measurement  noise,  (3)  actuator  errors,  (U)  system-to-system  parameter 
variations,  and  (5)  unpredictable  plant  disturbances. 

A design  approach  based  upon  separation  of  the  deterministic  and  stochas- 
tic aspects  of  the  control  problem  offers  a potential  means  for  meeting 
increased  Naval  requirements  associated  with  rapid  propulsion  system  response. 
This  design  procedure  involves  (l)  nonlinear  deterministic  feedback  control 
design,  (2)  nonlinear  estimator  design,  and  (3)  closed-loop  control  based  on 
feedback  of  estimated  system  state  variables  through  the  deterministic  control 
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logic.  The  analytical  techniques  required  for  development  of  nonlinear 
deterministic  control  logic  — in  particular,  piecewise-linear/piecewise- 
optimal  (PLPO)  procedures  — have  been  developed  by  UTRC  under  ONR  support 
(Ref.  5).  Procedures  for  nonlinear  estimator  design  using  nonlinear  filtering 
logic  with  PLPO-computed  gains  and  compensation  for  model-mismatch  have 
been  developed  and  evaluated  in  an  ONR-supported  study  at  UTRC  (Ref.  6). 
Therefore,  this  research  study  was  directed  toward  combining  the  determinis- 
tic control  and  stochastic  estimation  algorithms  and  evaluating  the  result- 
ing closed-loop  engine/estimator/controller  by  application  to  a nonlinear 
digital  dynamic  computer  simulation  of  the  FlOO/FUOl  turbofan  engine.  Mode- 
switching logic  was  to  be  developed  to  provide  smooth  transition  between  - 
small-signal  regulation  and  large-signal  transient  modes  of  estimator /con- 
troller operation.  Stability-in-the-large  as  well  as  stability-in-the-small 
was  to  be  determined.  Performance  of  the  nonlinear  FlOO/F^Ol  engine  simula- 
tion/estimator/ccntroller  was  to  be  evaluated  for  (l)  nominal  and  degraded 
engine  models,  ( 2 ) nominal  and  off-nominal  noise  statistics,  and  (3)  small- 
signal  as  well  as  large-signal  time  responses.  An  important  constraint  on 
the  closed-loop  estimator/controller  to  be  developed  will  be  that  its  compu- 
tational requirements  be  compatible  with  projected  airborne  digital  computer 
capabilities. 
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The  nonlinear  stochastic  control  problem  is  presented  in  detail  in 
Ref.  6.  Reference  6 sets  forth  (l)  the  general  nonlinear  stochastic  system 
structure,  (2)  the  control  objectives  and  constraints,  and  (3)  the  UTRC  design 
approach.  This  approach  is  based  upon  separating  the  stochastic  and  determin- 
istic aspects  of  the  control  problem.  The  stochastic  control  design  procedure 
involves  (l)  nonlinear  deterministic  feedback  control  design,  (2)  nonlinear 
estimator  design,  and  (3)  closed-loop  control  based  on  feedback  of  estimated 
engine  variables  through  the  deterministic  control  logic.  The  analytical 
techniques  for  defining  the  nonlinear  deterministic  control  (Ref.  5)  and 
nonlinear  estimator  (Ref.  6)  have  been  developed.  These  techniques  for  deter- 
ministic control  design  and  estimator  design  were  evaluated  by  application  to 
a nonlinear  digital  dynamic  FlOO/F^Ol  engine  simulation  (Refs.  3 and  6, 
respectively) . 

In  this  section  the  procedures  for  deterministic  control  and  stochastic 
estimator  design  are  reviewed.  Engine  dynamics  are  described  first.  Deter- 
ministic control  logic  is  discussed  next.  The  stochastic  estimation  algorithm 
is  then  set  forth. 
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Engine  Dynamics 


The  system  model  is  shown  in  Fig.  1 with  provision  for  the  estimation 
and  control  algorithms  included.  The  system  consists  of  a nonlinear  engine 
model  and  sensor  models.  Fuel  flow  driving  noise  is  included  to  account  for 
metering  valve  uncertainty.  Engine  state  variables  are  generated  through  the 
engine  dynamics  and  the  actual  engine  inputs.  These  state  variables  in  con- 
junction with  the  actual  inputs  govern  engine  output  response.  Engine  state 
variables  are  available  through  sensors  which  contain  inherent  lags.  These 
sensors  indicate  which  state  variables  are  measured.  Sensor  noise  is  included 
to  account  for  measurement  inaccuracies. 

A nonlinear  F100/FL01  turbofan  engine  simulation  (Ref.  7)  was  selected 
for  comparative  evaluation  of  the  control  and  estimation  algorithms.  The 
F100/FL01  afterburning  turbofan  is  a military  propulsion  system  incorporating 
the  latest  achievements  of  advanced  engine  technology  including  variable 
control  geometries  and  digital  electronic  supervisory  control.  The  variable 
geometries  as  well  as  the  nonlinear,  complex,  multivariable  interactions  among 
engine  input-output  variables  make  the  F100/FA01  engine  a particularly  chal- 
lenging application  of  modern  control  theory.  Dynamics  of  this  nonlinear 
simulation  are  represented  by  sixteen  state  variables.  A set  of  sixteen  engine 
state  variables  is  listed  in  Ref.  6. 
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The  complexity  of  a deterministic  control  and  stochastic  estimator 
designed  using  modern  control  and  estimation  theory  — in  particular,  linear 
quadratic  Gaussian  theory  — is  directly  related  to  the  order  of  the  linear 
engine  design  model  (i.e.,  the  number  of  states  used  in  the  engine  model 
employed  in  the  design  procedure).  Including  all  state  variables  in  the 
design  model  results  in  unnecessary  complexity  in  the  design  effort  as  well 
as  in  implementation  requirements  of  the  resulting  combined  estimator/con- 
troller.  It  is  therefore  desirable  to  reduce  the  order  of  the  full-state 
linearized  engine  model  by  retaining  only  significant  engine  dynamics  for 
estimator/controller  design. 

The  first  step  toward  determining  a reduced-order  linearized  engine 
model  for  control  and  estimator  design  is  to  select  the  engine  state  vari- 
ables. A fifth-order  design  model  was  chosen.  Engine  state  variables 
selected  are: 


• fan  turbine  inlet  temperature  (x^) 

• main  burner  pressure  (x^) 

• fan  speed  (x^) 


•compressor  speed  (x^) 

• afterburner  pressure  (x^). 


These  five  engine  states  are  available  through  sensors  modeled  by  first- 

order  lags  and  noise  sources.  Therefore,  in  addition  to  the  five  selected 

engine  states  there  are  five  sensed  states  — the  outputs  of  the  sensors. 

The  measurements  represent  noise-corrupted  sensed  engine  states.  The  sensed 

states  x,  ....  x and  roisy  measurements  z,  ....  zr  are: 

6 10  15 


sensed  x. 


z = noise-corrupted  x, 
1 o 


x„  = sensed  x 
7 2 


x0  = sensed  x„ 
0 3 


x„  = sensed  x, 

9 ^ 


x^  = sensed  xr 


noise-corrupted  x 
noise-corrupted  Xg 
noise-corrupted  x 
noise-corrupted  x^ 
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'he  engine  input  variables  chosen  are: 


• Jet  exhaust  area  (u^) 

• fan  inlet  guide  vanes  (u^) 

• compressor  variable  vanes  ( 

• main  burner  fuel  flow  (u^). 
Critical  engine  output  variables  are: 


) 


• fan  corrected  airflow  (y  ) 

• fan  stability  margin  (y^) 

• compressor  stability  margin  (y^) 

• thrust  (y^) 

• high  turbine  inlet  temperature  (y^). 

These  engine  output  variables  cannot  be  measured  directly  in  flight.  However 
the  time  evolution  of  these  critical  variables  determines  satisfactory  engine 
performance.  Accurate  in-flight  estimation  of  these  variables  will  lead  to 
more  effective  dynamic  control  (Ref.  8). 


Deterministic  control  and  stochastic  estimation  throughout  the  idle  to 
military  sea-level  static  operating  regime  (9  to  100  percent  thrust)  was 
investigated.  Table  I lists  values  of  the  steady-state  engine  variables  for 
sea-level  static  conditions  at  five  power  lever  angle  (PLA)  settings: 

PLA  = 20,  35,  ^7,  60  and  73  deg.  For  convenience,  all  engine  parameters 
except  fan  and  compressor  vane  positions  and  fan  and  compressor  stability  mar 
gins  have  been  normalized  to  1.0  at  PLA  = 73  deg.  The  vane  positions  are 
defined  as  ratios  of  their  maximum  positions  and  the  two  stability  margins 
are  given  as  ratios  of  one  (a  smaller  value  indicates  reduced  stability  mar- 
gin). 


The  next  step  in  the  control  and  estimator  design  procedure  is  to  define 
reduced-order  linearized  system  dynamics  at  a series  of  points  along  the 
steady-state  operating  line.  Several  methods  for  computing  reduced- order 
dynamics  are  described  in  Ref.  6.  Reduced-order  FlOO/F^Ol  engine  models 
were  calculated  by  identifying  reduced-order  dynamics  directly  from  input- 
output-state  data  (Refs.  5 and  9). 


Linearized  dynamics  at  a steady-state  operating  point  are  represented 


by 


““"M-  [-$&]  K£<” 

8y(t)  = C8x#(t)  + D(8u(t)  + C(t)) 


(1) 


8z(t)  = 


8xe(t) ' 
.SxjjitT 


+ Tj(t) 


where  the  vectors  xe , Xg , u,  y,  and  z represent  the  five  (5)  engine  states, 
five  (5)  sensor  states,  four  (4)  engine  inputs,  five  (5)  engine  outputs  and 
five  (5)  measurements,  respectively,  denotes  an  i x J null  matrix,  and 
1^  represents  a j order  identity  matrix.  The  vectors  £ (4xl)  and  n (5xl) 
represent  process  and  sensor  noise,  respectively.  The  constant  Afe  (5x5), 

Be  (5x4),  C 15x5)  and  D (5x4)  matrices  define  engine  perturbational  dynamics 
about  the  steady-state  operating  point.  Elements  of  the  matrix  AT  (5x5)  are 


(AT)jj  = Tj  ' * = J 

(Ar)jj  = Ot'^j 


12) 


where  t.  represents  the  sensor  time  constant  associated  with  the  i^b  engine 
state.  Sensor  time  constants  and  standard  deviations  of  sensor  and  fuel  flow 
driving  noise  are  presented  in  Table  II.  Variable  geometry  uncertainties 
were  assumed  negligible. 


Deterministic  Control 

A systematic  technique  for  deterministic  multivariable  nonlinear  system 
control  design  based  on  linear  quadratic  regulator  theory  — specifically, 
the  piecewise-linear/piecewise-optimal  (PLPO)  control  technique  — was 
developed  (Ref.  5).  The  deterministic  control  design  procedure  assumes  no 
uncertainties;  i.e.,  it  is  assumed  that  (l)  no  actuator  errors  exist,  (2) 
no  plant  disturbances  occur,  (3)  all  state  and  output  variables  are  measured 
perfectly,  and  (4)  actuator  and  plant  dynamics  and  parameters  are  known 
exactly.  Under  these  assumptions,  plant  state  and  output  variables  can  be 
determined  for  any  given  commanded  inputs. 

The  deterministic  PLPO  design  procedure  was  applied  to  the  PfcWA  digital 
computer  dynamic  simulation  of  the  F100/F401  engine.  A nonlinear  multi- 
variable  feedback  controller  was  defined  for  idle  to  military  sea-level 
static  engine  operation  (9  to  100  percent  thrust).  The  analytical  design 
involved  (l)  linearizing  the  F100/F401  engine  dynamics  about  the  five 
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steady-state  operating  points  between  idle  and  military  thrust,  (2)  applying 
linear  optimal  control  synthesis  methods  at  each  point,  and  (3)  combining 
the  five  optimal  linear  controllers  into  a single  nonlinear  controller  which 
had  feedback  gains  scheduled  with  high  compressor  speed. 

The  deterministic  control  dynamics  are  given  by 

xe(t)  t 

u(t)=  j L(N2)  dxe(T)  + /H(N2)(v(T)-Vc(T))dT  + U(0)  (3) 

*e  1°)  o 


where  the  L (^x5)  and  H (UxM  matrices  represent  proportional  plus  integral 

feedback  controller  gains,  respectively.  The  high  compressor  speed  is 

denoted  by  N_.  The  vectors  v and  v are: 

2 c 

v = fan  inlet  guide  vanes  v = commanded  v 

1 cl  1 


v 


2 


compressor  variable  vanes 


v _ * commanded  v„ 
c2  2 


v 


3 


thrust 


v „ = commanded  v„ 

c3  3 


v^  = high 


turbine  inlet  temperature  v ^ = commanded  v 


V 


The  controller  gain  matrices,  L and  K,  at  the  steady-state  operating  conditions 
(FLA  = 20,  35,  ^*7,  60  and  73  deg)  are  listed  in  Ref.  3.  Note  from  v and  v^ 
above  that  this  control  mode  requires  feedback  of  thrust  and  high  turbine 
inlet  temperature  — two  of  the  five  unmeasurable  output  variables  — in 
addition  to  the  engine  states.  Also,  engine  states  are  available  only 
through  noise-corrupted  sensors  with  inherent  lags.  However,  in  the  deter- 
ministic control  study,  engine  states,  thrust,  and  high  turbine  inlet  tem- 
perature were  assumed  to  be  known  exactly. 


Results  obtained  in  the  deterministic  dynamic  control  study  indicate 
the  ability  of  multivariable  feedback  control  logic  — using  feedback  of 
actual,  not  estimated,  variables  — to  provide  significantly  more  rapid 
propulsion  system  dynamic  response  than  that  associated  with  conventional 
control  logic  (Refs.  2 and  3). 


Estimation  Algorithm 

State  and  output  estimation  procedures  for  stochastic  nonlinear 
dynamic  systems  were  developed  and  evaluated.  The  stochastic  aspects  of  the 
problem  are  reintroduced  for  this  estimation  portion  of  control  system  design. 


G 


In  addition  to  process,  sensor,  and  model  uncertainties,  the  fact  that  all 
engine  variables  cannot  be  measured  and  that  any  measurement  is  subject  to 
sensor  errors  was  taken  into  account.  Model-mismatch  compensation  techniques 
were  employed  to  account  for  potential  mismatch  between  the  system  model  in 
the  filtering  algorithm  and  the  actual  nonlinear  system  model.  The  filtering 
algorithm  was  evaluated  by  application  to  the  nonlinear  digital  dynamic  com- 
puter simulation  of  the  FlOO/F^Ol  turbofan  engine.  An  important  constraint 
on  the  estimation  algorithm  is  that  its  computational  requirements  be  com- 
patible with  projected  airborne  digital  computer  capabilities. 

A large-signal  estimation  algorithm  based  on  Kalman  estimation  theory 
with  mo del- mis match  compensation  was  developed.  Results  obtained  in  a 
previous  UTRC  study  directed  toward  stochastic  small-signal  regulation  of 
nonlinear  multivariable  dynamic  systems  indicate  that  to  achieve  accurate 
estimation  Kalman  filtering  methodology  with  model-mismatch  compensation 
should  be  employed  (Ref.  8).  In  addition,  the  Ref.  8 study  showed  that 
improved  estimation  leads  to  improved  stochastic  regulation. 


Dynamics  of  nonlinear  filtering  logic  are^given  by 

,f|(*I(t),u(t),0)  + Kf|  (N2)(z(t)-  $,(t))  Kc(A2)  (z(r)  - fis(r))dr 
= + Kfz(N2)(z(t)-5#(t)) 

y(t)=  g(i*  (t),u(t),o) 


(*) 


where  (A)  denotes  the  estimate  of  the  variable  in  parentheses.  The  xe+  nota- 
tion denotes  full-state  representation  whereas  xe  denotes  reduced-state  repre- 
sentation. The  vector  function  f represents  reduced-order  engine  state 
dynamics.  The  Kalman  filter  matrices,  Kf^  (5x5)  and  Kf^  (5x5),  and  the  com- 
pensator gains,  Kc  (5x5),  are  calculated  off-line  at  the  five  (5)  operating 
points  based  on  the  fifth-order  linear  representations  of  the  nonlinear 
engine  simulation.  These  gains  are  computed  using  Kalman  filtering  theory 
as  well  as  model-mismatch  compensation  logic  (Ref.  6).  The  filter  gains  are 
tabulated  in  Ref.  6. 

Results  obtained  show  that  precise  estimation  was  achieved  using  the 
nonlinear  filtering  logic  with  model -mismatch  compensation.  Accurate  esti- 
mation of  key  engine  states  and  outputs  from  (1)  nominal-engine  data,  (2) 
degraded-engine  data,  and  (3)  engine  data  with  off-nominal  noise  statistics 
(that  is,  noise  statistics  different  from  the  statistics  assumed  in  calculat- 
ing the  filter  gains)  was  obtained  with  the  nonlinear  filtering  algorithm. 

In  addition,  the  computational  requirements  of  the  nonlinear  filtering  logic 
are  significantly  less  than  the  computational  requirements  of  an  extended 
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Kalman  filter.  The  major  reduction  in  computation  is  due  to  gain 
calculation  (Ref.  6).  Gain  calculation  for  an  extended  Kalman  filter 
designed  using  the  reduced-order  engine  model  and  not  including  on-line 
linearization  would  require  5325  multiplications  per  sampling  interval  and 
1» 865  additions  per  sampling  interval.  These  computations  are  based  on  ten 
states  (five  engine  and  five  sensed  states),  four  inputs,  and  five  measure- 
ments. On  the  other  hand,  gain  calculation  for  the  nonlinear  filtering 
algorithm  for  the  sea-level  static  flight  condition  considered  here  requires 
only  75  multiplications  per  sampling  interval  and  150  additions  per  sampling 
interval.  Including  sensor  model  mismatch  would  increase  the  required  multi- 
plications and  additions  by  25  and  50,  respectively.  These  results  imply 
that  the  nonlinear  estimation  logic  described  here  provides  a potentially 
viable  approach  to  estimation  of  key  engine  variables  under  realistic  operat- 
ing conditions. 
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COMBINED  ESTIMATION  AND  CONTROL  DESIGN 


In  this  section  combining  the  deterministic  control  and  stochastic 
estimator  into  a unified  feedback  controller  is  presented.  Closed-loop 
linearized  system  dynamics  are  introduced  first.  The  relationship  between 
control  and  estimator  dynamics  is  discussed.  Assessing  closed-loop  system 
stability  is  then  set  forth.  Mode-switching  logic  to  achieve  rapid  small- 
signal  as  well  as  large-signal  response  is  described  next.  Also,  mode- 
switching logic  to  achieve  rapid,  stable  degraded  engine  as  well  as  nominal 
engine  response  is  discussed.  In  the  final  part  of  this  section  com- 
putational requirements  for  the  combined  estimation  and  control  algorithm 
are  delineated. 


Closed-Loop  System  Dynamics 


Closed-loop  linearized  engine,  sensor,  filter,  and  control  dynamics  are 
given  by  the  differential  equations 
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where  w represents  linearized  engine  model  uncertainty. 


An  alternate  state-space  representation  of  the  closed-loop  system 
dynamics  is  obtained  by  defining  engine  state  and  sensor  estimation  errors 
(e-^  and  eg),  respectively,  by 

A 


«l  * *«"*• 
*2  * . 


(6) 


Substituting  x 


x - e. , x ■ x 
0 0 I s S 

closed-loop  linearized  system  dynamics 


e from  Eq.  (6)  into  Eq.  (5)  leads  to 


(5) 
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From  Eq.  (7)  it  can  be  seen  that  the  eigenvalues  of  the  closed-loop  system 
are  given  by  the  eigenvalues  of  the  matrices 


F„(9x9)A 


Ae  Be 

LAe  + HCe  LBe+HDe 
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-(AT-FKf2)  0 
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(10) 


That  is,  the  closed-loop  linearized  system  eigenvalues  are  (l)  the  eigen- 
values associated  with  the  deterministic  control  (Eq.  (8)),  (2)  the  eigen- 
values associated  with  the  sensor  dynamics  (Eq.  (9)),  and  (3)  the  eigen- 
values associated  with  the  stochastic  estimator  (Eq.  (10 ) ) . The  deter- 
ministic control  eigenvalues  can  therefore  be  adjusted  by  the  value  of  the 
feedback  control  gains  H and  L (which  depend  on  the  weighting  matrices 
within  the  control  performance  criterion)  while  the  estimator  eigenvalues 
can  be  adjusted  by  the  value  of  the  filter  gains  (which  depend  on  the  process, 
sensor,  and  model-mismatch  noise  intensity  matrices). 

If  the  eigenvalues  of  the  deterministic  control  are  dominant  (i.e., 
closer  to  the  Ju>  axis  in  the  left  half  plane)  the  transient  response  of  the 
closed-loop  system  in  the  stochastic  environment  will  be  similar  to  that  of 
the  deterministic  response.  However,  the  system  response  will  be  noise- 
corrupted  because  the  filter  will  pass  much  of  the  measurement  noise.  This 
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noise  acts  as  a disturbance  on  the  control  system.  Designing  the  filter/ 
control  so  that  the  control  eigenvalues  are  dominant  appears  to  be  applicable 
only  when  measurements  are  very  accurate.  On  the  other  hand,  if  the  eigen- 
values of  the  estimator  are  dominant,  closed-loop  system  transient  response 
will  be  different  from  the  deterministic  response.  Therefore,  the  controller 
should  be  designed  so  that  the  estimator  and  deterministic  control  eigen- 
values are  similar  if  the  system  measurements  are  noise-corrupted. 

Filter  and  control  gains  at  five  design  operating  points  (PLA  = 20,  35  > 
1*7,  60  and  73  deg)  were  calculated  previously.  Model-mismatch  uncertainty 
was  selected  in  the  estimator  design  so  that  neither  the  control  nor  filter 
eigenvalues  were  dominant.  Error  indices  — which  measure  the  normalized 
root  mean-square  error  between  linear  and  nonlinear  model  time  responses  — 
were  employed  as  a measure  of  the  linearized  model  uncertainties.  Table  III 
lists  the  eigenvalues  associated  with  the  sensor  dynamics,  with  the  deter- 
ministic control,  and  with  the  stochastic  estimator  at  PLA  = 73  deg.  This 
table  indicates  that  the  magnitudes  of  the  eigenvalues  associated  with  the 
stochastic  estimator  are  similar  to  the  magnitudes  of  the  eigenvalues  assoc- 
iated with  the  deterministic  control,  i.e.,  neither  the  estimator  nor  con- 
troller dynamics  are  dominant.  Similar  results  were  obtained  at  the  other 
design  points. 


Stability  Analysis 

Stability-in- the- small  at  design  points  is  assured  under  appropriate 
assumptions  for  filter  and  control  design  using  PLPO  techniques.  Closed- 
loop  system  stability-in-the-small  at  off-design  points  is  not  guaranteed 
due  to  the  interpolation  of  estimator  and  control  gains.  Stability  at  the 
off-design  points  was  determined  directly  from  computation  of  closed-loop 
system  eigenvalues  at  those  points.  To  assess  stability-in-the-small  at 
off-design  operating  points  (l)  open-loop  linearized  F100/FU01  engine  models 
were  identified  at  the  off-design  operating  points,  and  (2)  closed-loop 
stability-in-the-small  was  determined  from  the  linear  models  with  appropriate 
control  and  estimator  gains. 

To  identify  op>en-loop  linearized  models  for  determining  stability-in- 
the-small  at  off-design  operating  points  deterministic  input-output-state 
data  were  generated  by  the  nonlinear  FlOO/FUOl  engine  simulation  at  four 
off-design  steady-state  operating  conditions.  The  four  off-design  opera- 
ting points  selected  were  midway  as  a function  of  PLA  between  the  design 
operating  points  (that  is,  PLA  * 28,  Ul,  5U,  67  deg).  To  generate  the  data 
the  engine  inputs  were  sequentially  stepped  at  the  steady-state  off-design 
operating  conditions.  The  step  input  duration  was  one  second  with  one-half 
second  between  the  input  steps.  The  step  magnitudes  were  approximately 


jrrn  *?."•>  . -V* 


#7 


14 


(l)  five  percent  military  operating  point  on  Jet  exhaust  area,  (2)  five  and 
ten  percent  of  maximum  position  on  fan  inlet  guide  vanes  and  rear  compressor 
variable  vanes,  respectively,  and  (3)  five  percent  of  operating  point  on 
main  burner  fuel  flow.  These  input  perturbations  are  the  same  as  those 
perturbations  used  to  identify  linearized  models  at  the  design  operating 
points . 

Open-loop  linearized  models  for  determining  stability-in-the-small  were 
identified  at  the  four  off-design  operating  conditions  from  deterministic 
input-output-state  data  generated  by  the  nonlinear  FlOO/F^Ol  engine  simula- 
tion. The  least-squares  identification  technique  was  used  to  compute 
parameters  of  the  linearized  models  from  this  engine  data.  The  combination 
of  sequential  step  inputs  and  least-squares  parameter  identification  has  been 
an  effective  approach  for  defining  credible  linear  system  models  from  noise- 
free  data  when  all  variables  are  available. 

State  and  output  error  indices  at  the  four  off-design  operating  condi- 
tions are  listed  in  Table  IV.  These  error  indices  are  normalized  root-mean- 
square  differences  between  actual  and  calculated  engine  model  time  responses. 
Previous  UTRC  studies  indicate  that  identified  engine  dynamics  with  error 
indices  less  than  0.5  adequately  predict  small-signal  nonlinear  engine  per- 
formance for  control  and  estimator  design.  Table  TV  indicates  that  error 
indices  associated  with  the  identified  linear  models  at  the  off-design 
operating  points  are  much  smaller  than  0.5*  Actual  and  identified  engine 
model  state  responses  with  the  largest  and  smallest  error  indices  at  PLA  = 67 
deg  are  shown  in  Fig.  2.  Actual  and  identified  engine  model  output  responses 
with  the  largest  and  smallest  error  indices  at  PLA  * 67  deg  are  compared  in 
Fig.  3.  The  comparison  of  engine  responses  at  PLA  = 67  deg  is  representative 
of  results  at  the  other  three  off-design  operating  conditions.  These  time 
responses  in  conjunction  with  Table  IV  show  that  the  identified  linear  models 
accurately  predict  small-signal  nonlinear  engine  response. 

Closed-loop  stability-in-the-small  at  the  off-design  operating  con- 
ditions was  determined  from  the  identified  linear  models  with  appropriate 
control  and  estimator  gains.  Control  and  estimator  feedback  gain  matrices 
at  the  off-design  operating  points  were  determined  from  control  and  esti- 
mator gains  at  adjacent  design  operating  points  using  linear  interpolation 
with  compressor  speed  as  the  independent  variable.  (Linear  interpolation 
is  employed  to  schedule  filter  and  control  gains  as  a function  of  compressor 
speed  within  the  stochastic  control  logic  interfaced  with  the  nonlinear  F100/ 
FUOl  engine  simulation.)  Stability-in-the-small  at  off-design  operating 
points  was  demonstrated  by  computing  closed- loop  system  eigenvalues  (Eq.  (7)). 
Closed-loop  eigenvalues  at  PLA  * 67  deg  are  shown  in  Table  V.  This  table 
indicates  the  desirable  situation  where  neither  estimator  nor  controller 
dynamics  are  dominant.  Similar  results  were  obtained  at  the  other  three 
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off-design  operating  conditions.  In  addition,  magnitudes  of  the  eigenvalues 
at  off-design  operating  points  are  of  the  same  order  of  magnitude  as  the 
magnitudes  of  the  eigenvalues  at  the  adjacent  design  operating  conditions 
(see  Tables  III  and  V). 

The  control  and  estimator  synthesis  procedure  using  the  separation 
approach  does  not  guarantee  stability-in-the-large  for  nonlinear  system 
dynamics.  Stability-in-the-large  was  assessed  through  computer  simulation. 
Combined  engine  simulation /sensor  models/control  logic/  estimation  algorithm 
performance  as  well  as  stability-in-the-large  is  described  in  the  performance 
evaluation  section. 


Mode-Switching  Logic 

Mode-switching  logic  was  developed  to  provide  rapid  small-signal  thrust 
response.  The  deterministic  controller  feedback  gains,  optimized  for  large- 
signal  response,  are,  in  essence,  "tuned"  to  large  initial  errors  between 
commanded  and  actual  engine  variables.  When  these  errors  are  small  — as  in 
small-signal  transient  response  — these  deterministic  control  gains  are  too 
small  to  provide  rapid  small-signal  response.  To  obtain  rapid  small-signal 
as  well  as  rapid  large-signal  dynamic  response  mode-switching  logic  was 
developed.  The  mode-switching  logic  amplifies  the  commanded  engine  variable 
errors  (differences  between  estimated  and  commanded  engine  variables  which 
drive  the  deterministic  control  integral  path).  These  errors  are  amplified 
when  the  thrust  error  is  small  and  the  high  turbine  inlet  temperature  is  not 
near  its  maximum  limit.  The  amplification  gain  which  multiplies  the  commanded 
engine  variable  errors  is  given  by 


k = k|«2  (n) 

where  is  a function  of  the  thrust  error  and  Kj  is  a function  of  the 
estimated  high  turbine  inlet  temperature.  A schematic  diagram  of  the  mode- 
switching logic  to  achieve  rapid  small-signal  response  is  depicted  in  Fig.  i*. 
Figure  U shows  that  for  estimated  high  turbine  inlet  temperatures  less  than 
50  percent  military  operating  condition  Kg  is  equal  to  1.0  and  the  amplifica- 
tion gain  is  a function  of  the  thrust  error  only  (i.e.,  K = %).  For  esti- 
mated high  turbine  inlet  temperatures  greater  than  or  equal  to  75  percent 
military  operating  condition  (military  operating  condition  represents  the 
maximum  temperature  operating  limit).  Kg  is  equal  to  the  inverse  of  K^  and  K 
is  equal  to  1.0.  That  is,  when  the  high  turbine  inlet  temperature  is  at  or 
above  75  percent  of  its  safe  steady-state  operating  condition  the  errors  are 
not  amplified.  For  temperatures  between  50  and  75  percent  K2  varies  smoothly 
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between  1.0  and  the  inverse  of  and  K varies  smoothly  between  and  1.0. 

The  gain  K^  was  defined  as  a function  of  the  thrust  error,  as  shown  in  Fig. 

1*,  to  insure  smooth  transition  between  small-signal  and  large-signal  modes. 

If  the  thrust  error  is  greater  than  10  percent  military  thrust,  then  K = 1.0, 

Kg  =1.0  and  the  errors  to  the  controller  are  not  amplified.  For  thrust 
errors  ranging  from  10  percent  to  5 percent  varies  linearly  from  1.0  to 
25.0.  When  the  thrust  error  is  less  than  5 percent,  = 25.0.  The  slopes 
of  and  Kg  as  well  as  the  maximum  value  for  K^  were  chosen  to  yield  stable, 

rapid  small-signal  engine  thrust  response. 

In  addition,  logic  was  developed  to  achieve  stable,  rapid  degraded  engine 
performance.  The  deterministic  control  logic  steady-state  schedules  (requested 
steady-state  thrust,  high  turbine  inlet  temperature,  fan  and  compressor  vane 
positions  as  a function  of  PLA)  are  defined  for  the  nominal  engine.  As  the 
engine  degrades  these  schedules  can  no  longer  be  achieved.  The  critical 
variables  in  steady-state  operation  are  thrust  and  high  turbine  inlet  tempera- 
ture. A degraded  engine  operating  at  nominal  temperature  generally  will 
achieve  a thrust  less  than  nominal.  In  previous  UTRC  F100  deterministic  con- 
trol studies  the  deterministic  controller  provided  increased  thrust  by 
decreasing  area  to  its  lower  position  limit  thereby  saturating  the  area  con- 
trol. When  Jet  exhaust  area  is  unable  to  recover  from  saturation  unstable 
engine  operation  results.  To  prevent  area  saturation  as  steady-state  degraded 
engine  operation  is  reached  jet  exhaust  area  is  set  to  its  steady-state  value 
(as  a function  of  PLA)  as  engine  operation  reaches  steady-state.  When  the 
absolute  value  of  the  thrust  error  is  less  than  20  percent  military,  the 
normalized  area  is  ramped  at  0.5  sec  to  its  steady-state  value.  That  is. 

Jet  exhaust  area  is  calculated  by 


A j (t  + At)  = Aj(t)  + 0.5  At  for  'thrust  error|  < 0.2. 


(12) 


When  the  absolute  value  of  the  thrust  error  is  greater  than  or  equal  to  20 
percent  military  thrust  then  Jet  exhaust  area  is  computed  according  to 
Eq.  (3). 


Computational  Requirements 

To  code  the  deterministic  control  logic  on  a digital  computer  the  control 
equation  (Eq.  (3))  with  mode-switching  logic  (Eq.  (11)  and  Fig.  1*)  may  be 
represented  by 
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u (t  + At)  = u(t)  + L (N2  (t ))  A xe(t)  ♦Hc(N2(t))K(v(t)-vc(t)) 


(13) 


l 


I 

I 
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where  At  represents  the  known  sampling  rate,  and 

A xe(t)  = x^t)-  xe(t-At) 


Hc  = hA1  . 


(lM 


Linear  univariate  interpolation  for  determining  control  gains  L and  Hc  between 
operating  points  based  on  Np  as  well  as  to  compute  the  mode-switching  logic 
gains  and  as  functions  of  thrust  error  and  high  turbine  inlet  tempera- 
ture, respectively,  is  employed  in  the  estimator/controller.  That  is. 


9ke=nrMN2-(N2U  +(9ke)| 


(^ke^  i + i “ (^ke)i 
' (n2)  i + i " (N2)| 

where  i represents  a steady-state  operating  condition  and  ( II  ) . < 11  < ( K ) 

2 l — 2 — 2 i+p 

Control  computational  requirements  (from  Eqs.  (13)  and  (15))  consist  of 
multiplications  and  additions.  To  calculate  the  control  and  mode-switching 
logic  gains  (Eq.  (15))  38  multiplications  and  76  additions  per  sampling 
interval  are  required.  To  calculate  updated  inputs  (Eq.  (13))  ^1  multiplica- 
tions and  l»5  additions  per  sampling  interval  are  needed.  Therefore,  the 
deterministic  control  with  mode-switching  logic  requires  79  multipli- 
cations and  121  additions  each  sampling  interval  to  compute  the  engine 
inputs. 

Computational  requirements  for  the  nonlinear  filtering  logic  are  dis- 
cussed in  detail  in  Ref.  6.  Filter  computational  requirements  depend  upon 
(l)  filter  gain  calculations  and  (2)  the  engine  model  employed  in  the  filter. 
The  gain  calculations  require  only  75  multiplications  and  150  additions  per 
sampling  interval.  These  gain  calculations  are  significantly  less  than  those 
required  by  the  conventional  extended  Kalman  filter  (5185  multiplications  and 
hjk5  additions  per  sampling  interval.)  A minimum  complexity  nonlinear  engine 
model  should  be  selected  for  the  filtering  algorithm  to  reduce  model  computa- 
tional requirements.  Results  indicate  that  model -mismatch  compensation  logic 
is  able  to  compensate  for  differences  between  the  engine  model  employed  in  the 
filter  and  the  actual  engine.  The  filter  computations  (dependent  on  model  used 
in  filter)  added  to  the  control  computations  represent  the  computations 
required  to  implement  the  nonlinear  stochastic  control  algorithm. 
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In  this  section  stochastic  control  performance  is  evaluated.  To  evaluate 
closed-loop  engine  performance  including  stability-in-the-large  the  nonlinear 
FIOD/'FLQI  engine  simulation  was  interfaced  with  the  nonlinear  stochastic  feed- 
back control  algorithm.  The  closed-loop  engine/stochastic  control  structure 
is  shown  in  Fig.  1.  The  overall  closed-loop  system  consists  of  (l)  the  non- 
linear FlOC/F^Ol  engine  simulation,  (2)  sensor  models,  (3)  deterministic  PLPC 
control  logic  with  node-switching,  arid  (k)  nonlinear  filtering  logic  with  model- 
mismatch  compensation.  The  stochastic  controller  performance  when  interfaced 
with  a nominal  engine  is  discussed  first.  Small-signal  as  well  as  large- 
signal  response  is  evaluated.  Degraded-engine  performance  is  then  described. 
Engine  performance  when  measurement  and  fuel  flow  noise  statistics  are  differ- 
ent from  these  statistics  assumed  in  estimator  design  is  set  forth  in  the 
final  part  of  this  section. 

Nominal  Engine 

Stochastic  engine  accelerations  from  idle  (PLA  = 20  deg)  to  military 
(PLA  = 73  deg)  thrust  levels  were  computed  using  the  closed-loop  nominal 
nonlinear  engine,  sensors  with  the  nominal  noise  statistics,  filter,  and  control. 
These  large-signal  accelerations  were  compared  with  the  deterministic 
accelerations  (nominal  engine/controller  with  no  noise  and  no  estimator)  of 
the  Ref.  3 study.  Results  indicate  that  deterministic  and  stochastic  time 
responses  are  almost  identical.  Stochastic  thrust  response  is  compared  with 
the  deterministic  thrust  response  in  Fig.  5-  Figure  5 illustrates  that 
stochastic  thrust  response  for  the  nominal  engine  with  nominal  noise  statis- 
tics is  as  rapid  as  deterministic  thrust  response.  As  desired,  there  is  no 
difference  in  nominal -engine  stochastic  and  deterministic  thrust  responses  due 
to  the  stochastic  controller  mode-switching  logic  to  accommodate  engine 
degradation.  This  mode-switching  logic  was  not  required  in  the  deterministic 
control  since  nc  uncertainties,  including  engine  degradation,  were  considered 
in  the  deterministic  design  procedure.  The  stochastic  engine  response  of  the 
other  engine  variables  was  very  similar  to  the  deterministic  response.  The 
similar  stochastic  and  deterministic  engine  performance  results  from  the 
accurate  estimation  achieved  by  the  nonlinear  filtering  logic.  The  mean  and 
standard  deviation  of  state  and  output  estimation  errors  (where  the  estimation 
error  = actual -estimated  response)  are  shown  in  Table  VI.  Table  VI  indicates 
that  the  average  error  is  zero  to  at  least  three  significant  figures  for  all 
state  and  output  estimates.  The  maximum  standard  deviation  is  0.002  in  fan 
turbine  inlet  temperature. 

To  evaulate  closed-loop  small-signal  system  performance  stochastic 
engine  response  to  a PLA  step  from  20  to  23  deg  was  generated. 
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Small-signal  actual  engine  model  thrust  response  using  control  logic  vith  and 
without  mode-switching  logic  is  compared  in  Fig.  6.  Figure  6 indicates  that 
the  small-signal  thrust  response  with  mode-switching  logic  incorporated  into 
the  PLPO  controller  is  significantly  faster  than  the  thrust  response  without 
mode-switching  logic.  For  example,  the  11.9  percent  thrust  point  (90  percent 
of  the  requested  change  in  thrust)  is  reached  in  3.0  sec  using  mode-switching 
logic  versus  h.2  sec  when  mode-switching  logic  is  not  employed.  Similar  rapid 
small-signal  thrust  response  was  achieved  using  the  stochastic  controller 
with  mode-switching  logic  at  other  off-design  as  well  as  design  points. 

Results  obtained  indicate  that  (l)  the  controller  requires  mode-switching 
logic  to  obtain  rapid  small-signal  as  well  as  rapid  large-signal  dynamic 
response  and  (2)  mode-switching  logic  provides  smooth  transition  between  the 
large-  and  small-signal  modes  of  operation. 

Degraded  Engine 

To  evaluate  mode-switching  control  logic  developed  to  accommodate  engine 
degradation,  degraded-engine  response  to  a step  in  PLA  from  20  to  73  deg  was 
generated  using  the  stochastic  controller  without  and  then  with  the  mode- 
switching  control  logic  (Eq.  (12)).  Engine  degradation  was  accomplished  by 
decreasing  fan  and  compressor  efficiencies  four  (b)  percent  and  one  (l) 
percent,  respectively.  Large-signal  engine  performance  to  a step  in  PLA  from 
20  to  73  deg  was  not  satisfactory  using  the  controller  without  mode-switching 
logic.  The  thrust  did  not  reach  military  operation,  the  fan  surged,  and  the 
turbine  overtemperatured.  This  poor  engine  performance  results  from  the  fact 
that  steady-state  schedules  (requested  steady-state  thrust,  high  turbine 
inlet  temperature,  fan  and  compressor  vane  positions  as  a function  of  PLA) 
are  defined  for  the  nominal  engine.  As  the  engine  degrades,  these  schedules 
can  no  longer  be  achieved. 
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Degraded-engine  performance  in  conjunction  with  key  control  analysis 
results  (Ref.  10)  were  analyzed.  Engine  response  showed  that  within  2.5  sec 
vane  positions  and  high  turbine  inlet  temperature  were  at  their  requested 
values.  However,  the  thrust  level  was  ninety-four  percent  military  thrust 
(requested  value  for  thrust  at  PLA  = 73  deg  is  100  percent  military  thrust). 
Analysis  of  engine  inputs  showed  that  the  control  logic  was  driving  jet 
exhaust  area  to  its  minimum  position  limit  to  drive  the  thrust  error 
( estimated-requested  thrust)  to  zero.  Key  control  analysis  indicates  that 
Jet  exhaust  area  is  a key  control  input  for  steady-state  as  well  as  transient 
thrust  response. 

Degraded-engine  performance  with  the  mode-switching  logic  to  accommodate 
engine  degradation  was  satisfactory.  Stochastic  degraded  and  nominal  engine 
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thrust  responses  are  compared  in  Fig.  7.  Figure  7 shows  that  degraded  engine 
response  is  only  slightly  slower  than  the  nominal  engine  response  until  90 
percent  military  thrust  is  reached  (2  sec).  The  nominal  engine  reaches  100 
percent  military  thrust  at  2.5  sec.  On  the  other  hand,  the  degraded  engine 
does  not  reach  100  percent  military  thrust  until  5 sec . 

The  stochastic  degraded  engine  response  is  compared  with  deterministic 
degraded  engine  response  in  Fig.  8 (deterministic  response  is  the  response 
achieved  by  feedback  of  actual  not  estimated  variables).  Figure  8 shows 
that  stochastic  and  deterministic  responses  are  very  similar.  Therefore, 
differences  in  nominal  and  degraded  engine  response  with  the  stochastic  con- 
troller (that  is,  the  relatively  slow  degraded-engine  thrust  response  from  90 
to  100  percent  military  thrust)  results  from  degraded  engine  performance 
limitations,  not  from  estimation  errors  in  engine  variables. 

Estimation  error  statistics  (mean  and  standard  deviation)  for  degraded 
engine  response  using  the  stochastic  controller  are  shown  in  Table  VII. 

Table  vn  indicates  that  very  accurate  estimates  of  engine  variables  Eire 
achieved  with  the  nonlinear  filtering  logic.  The  absolute  average  error  is 
0.002  for  fan  and  high  turbine  inlet  temperature  estimates,  0.001  for  main 
burner  pressure  and  fan  and  compressor  stability  margin  estimates,  and  zero 
to  three  significant  figures  for  all  other  state  and  output  estimates.  The 
maximum  standard  deviation  is  0.002  in  fan  turbine  inlet  temperature;  the 
standard  deviation  is  zero  to  three  significant  figures  in  fan  and  compressor 
speed  and  0.001  for  all  other  state  and  output  estimates. 

Evaluation  of  degraded  engine  performance  indicates  that  (l)  mode- 
switching logic  to  accommodate  engine  degradation  is  required  for  stable, 
rapid  dynamic  engine  response  and  (2)  model-mismatch  compensation  within  the 
nonlinear  filtering  algorithm  results  in  accurate  steady-state  as  well  as 
transient  estimation. 


Off-Nominal  Noise  Statistics 

The  effect  of  off-nominal  noise  statistics  on  stochastic  closed-loop 
engine  performance  was  assessed.  Noise-corrupted  engine  input  and  measure- 
ment data  with  noise  statistics  different  from  those  assumed  in  calculating 
filter  gains  were  generated  for  the  degraded  as  well  as  nominal  engine. 
Standard  deviations  of  the  measurement  noise  were  doubled;  whereas,  fuel  flow 
driving  noise  standard  deviation  was  halved.  Engine  response  with  off-nominal 
measurement  and  input  noise  was  the  same  as  engine  response  with  nominal 
measurement  and  input  noise  statistics.  (See  Fig.  5»  stochastic  nominal 
engine  and  Fig.  7,  stochastic  degraded  engine  for  thrust  response.)  These 
results  indicate  that  the  nonlinear  stochastic  control  is  insensitive  to  a 
mismatch  between  actual  and  nominal  statistics. 
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Results  obtained  in  this  study  show  that  rapid,  stable  engine  performance 
is  achieved  using  the  stochastic  controller  with  model-mismatch  compensation 
and  mode-switching  logic.  Satisfactory  thrust  response  was  achieved  with  the 
nonlinear  stochastic  controller  for  (l)  nominal  engine  with  nominal  noise 
statistics,  (2)  degraded  engine  with  nominal  noise  statistics,  (3)  nominal 
engine  with  off-nominal  noise  statistics,  and  (M  degraded  engine  with  off- 
nominal  noise  statistics.  In  addition,  the  controller  provided  rapid,  stable 
small-signal  as  well  as  large-signal  dynamic  response. 


RESULTS  AND  CONCLUSIONS 


1.  The  nonlinear  stochastic  feedback  control  algorithm  was  interfaced  with  a 
nonlinear  FlOO/Fl+Ol  engine  simulation.  The  overall  closed-loop  system 
consists  of  (l)  the  nonlinear  FlOO/F^Ol  engine  simulation,  (2)  sensor 
models,  (3)  a deterministic  piecewise-linear/piecewise-optimal  (PLPO) 
control  algorithm  with  mode-switching  logic,  and  (l*)  a nonlinear  esti- 
mator with  model-mismatch  compensation.  Mode-switching  logic  was  devel- 
oped to  provide  smooth  transition  between  small-signal  and  large-signal 
engine  operation.  Also,  control  logic  to  accommodate  engine  degradation 
was  defined.  Model-mismatch  compensation  was  employed  to  account  for 
mismatch  between  the  system  model  within  the  filtering  logic  and  the 
actual  system  model.  Engine  variables  required  for  feedback  through  the 
deterministic  control  logic  include  unmeasurable  variables  (thrust  and 
high  turbine  inlet  temperature)  as  well  as  measured  variables  (fan  and 
compressor  speeds,  main  burner  and  afterburner  pressures,  and  fan  turbine 
inlet  temperature).  The  nonlinear  stochastic  controller  (deterministic 
PLPO  control/nonlinear  estimator)  was  evaluated  throughout  the  idle  to 
military  sea-level  static  operating  regime  (9  to  100  percent  thrust). 

(a)  Closed-loop  engine/control/filter  stability-in-the-small  as  well  as 
stability-in-the-large  was  assessed  at  design  and  off-design  oper- 
ating points.  Stability-in-the-small  at  design  points  is  assured 
under  appropriate  assumptions  using  the  design  procedures  employed 
in  this  study.  Stability-in-the-smell  at  off-design  operating  con- 
ditions was  demonstrated  from  computation  of  closed-loop  system 
eigenvalues  at  selected  off-design  operating  conditions.  Stability- 
in-the-large  was  shown  through  computer  simulation. 

(b)  Small-signal  engine  performance  using  control  logic  with  and  without 
mode-switching  logic  was  compared.  Small-signal  thrust  response 
with  mode-switching  logic  incorporated  into  the  controller  was  sig- 
nificantly faster  than  thrust  response  without  mode-switching  logic. 
Also,  mode-switching  logic  provided  a smooth  transition  between 
small-  and  large-signal  operation. 

(c ) Nominal  stochastic  engine  performance  (nominal  engine  model/control/ 
estimator /nominal  fuel  flow  and  measurement  uncertainties)  was  compared 
with  nominal  deterministic  engine  performance  (nominal  engine  model/ 
control/no  estimator  or  system  uncertainties).  Stochastic  thrust 
response  was  as  rapid  as  deterministic  thrust  response  for  both 
large-  and  small-signal  operation.  Stochastic  engine  response  of 
other  critical  variables  was  very  similar  to  the  deterministic 
response.  This  similarity  in  deterministic  and  stochastic  perfor- 
mance was  achieved  by  assuring  that  neither  control  nor  estimator 
dynamics  were  dominant. 


23 


(d)  The  stochastic  control  algorithm  was  also  evaluated  by  application  to 
a degraded  F100/F*+01  engine  simulation.  Engine  degradation  was 
achieved  through  decreases  in  fan  and  compressor  efficiencies  repre- 
sentative of  realistic  operating  conditions.  Stochastic  degraded 
engine  response  was  slower  than  nominal  engine  response.  Comparison 
of  stochastic  and  deterministic  engine  performance  shows  that  the 
slower  degraded  engine  response  results  from  degraded  engine  perfor- 
mance limitations  rather  than  from  estimation  errors  in  engine 
variables . 

(e)  Performance  of  the  stochastic  control  algorithm  was  evaluated  for 
measurement  and  fuel  flow  metering  noise  statistics  different  from 
those  assumed  in  the  design  procedure.  Engine  response  with  off- 
nominal  noise  was  the  same  as  the  response  with  nominal  noise.  That 
is,  the  closed-loop  system  response  was  insensitive  to  a realistic 
mismatch  in  actual  and  assumed  noise  statistics. 

Results  obtained  indicate  that  the  nonlinear  stochastic  controller  design 
procedures  employed  herein  provide  a viable  approach  for  gas  turbine  con- 
trol design. 

(a)  Stability-in-the-small  and  in-the-large  was  established.  Rapid, 
stable  thrust  response  was  achieved  (l ) for  small-signal  as  well  as 
large-signal  operation,  (2)  for  a nominal  engine  model,  (3)  for  a 
degraded  engine  model,  and  (U)  for  off-nominal  noise  statistics. 

(b)  Computational  requirements  for  the  stochastic  controller  depend  upon 
(l)  control  gain  calculations,  (2)  control  computations  to  update 
engine  inputs,  (3)  filter  gain  calculations,  and  (**)  the  engine 
model  employed  in  the  filter.  The  control  and  filter  gain  calcula- 
tions as  well  as  the  input  update  calculations  require  only  15** 
multiplications  and  271  additions  per  sampling  interval.  A minimum 
complexity  engine  model  should  be  selected  for  the  filtering  algo- 
rithm to  reduce  model  computational  requirements.  Model-mismatch 
compensation  logic  should  be  employed  to  compensate  for  differences 
between  the  engine  model  in  the  filter  and  the  actual  engine. 


RECOMMENDATIONS  FOR  FUTURE  STUDY 
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Results  of  this  study  indicate  that  a multivariable  stochastic  controller 
designed  using  modern  state-space  estimation  and  control  techniques  is  a 
feasible  approach  to  control  of  complex,  nonlinear  variable  geometry  gas 
turbine  engines.  Studies  that  represent  a logical  broadening  and  analytical 
extension  of  the  research  reported  herein  are: 

• Exhaust  gas  re ingest ion:  Apply  the  control  and  estimation  methodology 
developed  in  this  study  to  define  and  evaluate  control  and  estimation 
logic  necessary  to  maintain  thrust  during  reingestion  of  engine  exhaust 
gases.  Problems  resulting  from  hot  gas  ingestion  assume  increasingly 
greater  importance  as  the  Navy  focuses  its  attention  upon  future  large- 
scale  deployment  of  high-speed  tactical  V/STOL  aircraft. 

• Hyb ri d eval uat ion:  Evaluate  the  nonlinear  stochastic  controller  on  a 
real-time  detailed  FlOO/FhOl  hybrid  simulation.  The  control  law  would 
be  coded  on  a digital  minicomputer.  Hardware  requirements,  system 
accuracy,  stability  as  a function  of  sampling  rate,  and  control  com- 
putation time  would  be  delineated.  This  evaluation  would  be  directed 
toward  eventual  engine  test  of  the  stochastic  controller. 

• Diagnostics : Develop  and  evaluate  diagnostic  algorithms  for  predicting 
and  identifying  engine  degradation  using  the  nonlinear  estimator.  The 
nonlinear  estimator  with  model-mismatch  compensation  achieves  accurate 
estimates  of  degraded  engine  variables.  Power  spectral  analysis  of 
the  error  between  estimated  variables  and  nominal  steady-state  values 
could  be  employed  to  identify  engine  degradation.  Results  of  the 
analyses  could  be  employed  for  maintenance  scheduling.  In  addition, 
adaptive  modulation  of  control  gains  based  on  the  emount  of  identified 
degradation  could  be  employed  to  achieve  rapid,  safe  acceleration  as 
well  as  accurate  steady-state  regulation  despite  engine  degradation. 

• Performance-seeking  logic:  UTRC  studies  carried  out  under  NSF  support 
(Ref.  M have  resulted  in  successful  development  of  adaptive  control 
logic  for  in-flight  minimization  of  thrust  specific  fuel  consumption 
(TSFC).  The  developed  performance-seeking  algorithm  was  evaluated  on 
the  FlOO/F^Ol  engine  simulation  (Ref.  7)  using  actual,  not  estimated 
variables.  The  nonlinear  estimator  described  herein  provides  a 
potential  means  for  estimating  unmeasurable  engine  variables  for  use 
within  the  performance-seeking  logic.  A program  should  be  conducted 

to  interface  the  estimator  with  the  developed  TSFC  in-flight  optimization 
algorithm.  TSFC  improvements  over  a range  of  representative  flight 
conditions  should  then  be  determined. 
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LIST  OF  SYMBOLS 


Constant  5x5  matrix  used  to  describe  linearized  engine  dynamics 

Jet  exhaust  area,  normalized 

5x5  matrix  used  to  represent  sensor  dynamics 

Element  in  the  ith  row  and  column  of  the  matrix  AT 

Constant  5x^  matrix  used  to  describe  linearized  engine  dynamics 

Constant  5x5  matrix  used  to  describe  linearized  engine  dynamics 

Constant  ^xh  matrix  used  to  describe  linearized  engine  dynamics 

5x1  error  vector  used  to  represent  the  difference  between  actual 
and  estimated  engine  state  perturbations 

5x1  error  vector  used  to  represent  the  difference  between  sensed 
and  estimated  sensor  state  perturbations 

Constant  9x9  matrix  used  to  represent  engine  control  dynamics 

Constant  5x5  matrix  used  to  represent  sensor  dynamics 

Constant  15x15  matrix  used  to  represent  stochastic  estimator 
dynamics 

Nonlinear  5x1  vector  function  — mathematical  representation  of 
rate-of- change  of  reduce.'  order  engine  state  vector 

Nonlinear  5x1  vector  function  — mathematical  representation  of 
engine  output  vector 

General  gain 

UxU  optimal  integral  feedback  regulator  gain  matrix 
HAt 

5x5  identity  matrix 

Amplification  gain  in  mode-switching  logic 

Component  of  mode-switching  logic  gain  — function  of  thrust  error 


LIST  OF  SYMBOLS  (Continued) 


Component  of  mode-switching  logic  gain  — function  of  high  turbine 
inlet  temperature 

5x5  compensator  gain  matrix  for  estimator  with  mo del -mis match 
compensation 

5x5  upper  partition  of  Kalman  filter  gain  matrix  for  estimator  with 
model-mismatch  compensation 

5x5  lower  partition  of  Kalman  filter  gain  matrix  for  estimator  with 
model-mismatch  compensation 

1x5  optimal  integral  feedback  regulator  gain  matrix 
Slope  associated  with  linear  gain  interpolation 
High  compressor  speed,  normalized 
i x,’  null  matrix 


Power  lever  angle,  degrees 
time,  sec 
Sample  time,  sec 
1x1  control  vector 

1x1  vector  used  to  represent  engine  variables  whose  steady-state 
values  are  specified 

1x1  vector  used  to  represent  the  commanded  values  of  v 
5x1  vector  used  to  represent  filter  model-mismatch 
10x1  system  state  vector 


5x1  engine  state  vector 
16x1  engine  state  vector 


5x1  sensor  state  vector 


5x1  engine  output  vector 


LIST  OF  SYMBOLS  (Continued) 


z 5x1  engine  measurement  vector 

5(  ) Perturbational  value  of  quantity  in  parenthesis 
n 5x1  sensor  noise  vector 

i l*xl  actuator  (process)  noise  vector 

Sensor  time  constant  associated  with  the  i^h  engine  state 
( ) Estimated  value  of  quantity  in  parenthesis 

(')  Time  derivative  of  quantity  in  parenthesis 

d(  ) Differential  of  quantity  in  parenthesis 

hi  ) Finite  increment  of  quantity  in  parenthesis 

h Equals  by  definition 
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TABLE  I 


NORMALIZED  STEADY- STATE  ENGINE  PARAMETERS 
AS  FUNCTIONS  CF  POWER  LEVER  ANGLE 


Sea-Level  Static  Operating  Condition 
Idle  Condition:  PIA  = 20  deg 

Military  Condition:  PIA  = 73  deg 


Engine  Parameters 

Power  Lever  Angle,  PLA-deg 

Type 

Parameter 

20 

35 

47 

60 

73 

Far.  Airflow 

KEI 

0.62 

0.73 

o.es 

1.0 

Fan  Stability  Margin 

Eel 

0.14 

0.18 

0.12 

0.12 

Output 

Compressor  Stability  Margin 

K&a 

0.18 

0.18 

0.19 

0.19 

Thrust 

0.09 

0.35 

0.52 

0.72 

1.0 

High  Turbine  Inlet  Temperature 

0.56 

0.73 

0.81 

0.90 

1.0 

Fan  Turbine  Inlet  Temperature 

0.57 

0.72 

0.81 

0.90 

1.0 

Main  Burner  Pressure 

0.21 

0.45 

0.59 

0.77 

1.0 

State 

Fan  Speed 

0.49 

0.74 

O.83 

0.91 

1.0 

Conpressor  Speed 

0.69 

0.83 

0.88 

0.93 

1.0 

Afterburner  Pressure 

0.40 

0.55 

0.66 



0.82 

Control 


Jet  Exhaust  Area 
Fan  Inlet  Guide  Vanes 
Conpressor  Variable  Vanes 
Main  Burner  Fuel  Flow 


O.98 

•0.50 

.1.11 

0.12 


0.98 

-0.50 

-0.39 

0.33 


O.98 

•0.50 

■0.17 

0.48 


0.98 

-0.34 

0.04 

0.70 


1.0 

-0.08 

0.20 

1.0 
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TABLE  III 


0 .297 


TABLE  IV 


ERROR  INDICES  ASSOCIATED  WITH  LINEAR 
FI00/F401  ENGINE  MODELS 


ERROR  INDEX  = 


I (ACTUAL  RESPONSE-IDENTIFIED  RESPONSE)2 


N , .2 

£ (ACTUAL  RESPONSE) 


N = Number  of  Data  Points 


Fan  Turbine  Inlet  Temperature 
Main  Burner  Pressure 
Fan  Speed 
Compressor  Speed 
Afterburner  Pressure 


Fan  Corrected  Airflow 
Fan  Stability  Margin 
Compressor  Stability  Margin 
Thrust 

High  Turbine  Inlet  Tanperature 


3 
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V. 


TABLE  VI 


fioo/fLoi  estimation  error  statistics  associated 

WITH  NOMINAL  ENGINE 


Estimation  Error  _A  Actual-Estimated  Response 


Variable 

Type 

Variable 

Average 

Error 

Standard 

Deviation 

of  Error 

Fan  Turbine  Inlet  Temperature 

0.000 

0.002 

Main  Burner  Pressure 

0.000 

0.001 

State 

Fan  Speed 

0.000 

0.000 

Compressor  Speed 

0.000 

0.000 

Afterburner  Pressure 

0.000 

0.001 

Fan  Corrected  Airflow 

0.000 

0.000 

Fan  Stability  Margin 

0.000 

0.001 

Output 

Compressor  Stability  Margin 

0.000 

0.001 

Thrust 

0.000 

0.001 

High  Turbine  Inlet  Temperature 

0.000 

0.001 

WITH  MODEL-MISMATCH 
COMPENSATION  LOGIC 


COMPRESSOR  SPEED  MAIN  BURNER  PRESSURE 


thrust  compressor  stability  margin 

(ERROR  INDEX  = 0.018)  , (ERROR  INDEX  = 0.290) 


NORMALIZED  F100/F401  ENGINE  MODEL  OUTPUT  RESPONSES  WITH 
LARGEST  AND  SMALLEST  ERROR  INDICES 

POWER  LEVER  ANGLE  - 67  DEG 


■ ACTUAL  NONLINEAR  MODEL  RESPONSE 
__  — IDENTIFIED  LINEAR  MODEL  RESPONSE 


TT-IO -63-2 


PERCENT  MILITARY 
OPERATING  CONDITION 


THRUST 


THRUST 


NORMALIZED  STOCHASTIC  NOMINAL  F100/F401  SMALL-SIGNAL  THRUST  RESPONSE 


STEP  IN  POWER  LEVER  ANGLE  FROM  20  TO  23  DEG 

ACTUAL  RESPONSE  WITH  MODE-SWITCHING  LOGIC 

— — — ACTUAL  RESPONSE  WITHOUT  MODE-SWITCHING  LOGIC 


TIME-SEC 


77-12-tO-l 
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